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In this paper we calculate some exact solutions of the grand partition functions for quantum 
gases in confined space, such as ideal gases in two- and three-dimensional boxes, in tubes, in annular 
containers, on the lateral surface of cylinders, and photon gases in three-dimensional boxes. Based 
on these exact solutions, which, of course, contain the complete information about the system, 
we discuss the geometry effect which is neglected in the calculation with the thermodynamic limit 
V ~* oo, and analyze the validity of the quantum statistical method which can be used to calculate 
the geometry effect on ideal quantum gases confined in two-dimensional irregular containers. We 
also calculate the grand partition function for phonon gases in confined space. Finally, we discuss 
the geometry effects in realistic systems. 
PACS numbers: 05.30.-d, 68.65.-k 

I. INTRODUCTION 

In statistical mechanics, for seeking the sum over all possible states, we always take the approximation that the 
volume of the system tends to infinity. In so doing we have lost all information about the geometry property of the 
system because in such an approximation the spectrum of single-particle states is continuous while the total number 
of states is independent of the shape of the boundary and simply proportional to the volume of the system. However, 
the properties of some systems found recently are shape dependent and sensitive to the topology [1-3]. In Ref. [4], we 
developed a method for calculating the effects of boundary and topology on ideal Bose and Fermi gases in confined 
space. 

In the thermodynamic limit, one can replace the summation over states by an integral: — > / dsp{e), where 



p(e) = y 



(1) 



is the density of states of a d-dimensional gas of dispersion E — p^/ {2m). However, if the system is enclosed in a finite 
volume, the structure of the phase space will be changed. In confined space, the spectrum of single-particle states 
will depend on the shape of the boundary. This is important if the thermal wavelength of particles is comparable to 
the size of the container, that is to say the replacement of the summation with the integral with the density of states 
Eq.(l) is a good approximation only when the volume is large enough so that the particle can not feel the boundary. 

For taking into account geometry effects, in the preceding paper [4] we proposed that for the case of a two- 
dimensional ideal quantum gas confined in an irregular container, the density of states, after neglecting the topological 
contribution, can be written as 

where S is the area and L the perimeter of the container. Moreover, for the case of a three-dimensional ideal quantum 
gas in a long tube, of which all transverse cross-sections keep the same, the density of states is 

^ „(27rm)3/2 1 27rm 
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S and L here denote the area and perimeter of the transverse cross-section of the tube, respectively, and the length 
of the tube is made sufficiently large so that the z-componcnt of the momentum pz can be considered to be 
continuous. The second terms of Eqs. (2) and (3) describe the contributions from the boundary. 

In the present paper, for analyzing geometry effects, we will calculate some exact solutions for ideal gases and 
photon gases in various kinds of confined space, and provide an approximate result for phonon gases in confined 
space. By comparing these exact and approximate solutions with the result obtained in the thermodynamic limit 
V DO, we can extract the information about boundary shape and topology. Furthermore, the validity of the method 
presented in Ref. [4] can be further justified by comparing it with these exact solutions. 

It can be found that the corrections to the standard results which are on the order X/L where A is the thermal 
wavelength and L denotes the linear size of the system. At the end of this paper, we will discuss such geometry effects 
in realistic systems. 

In Sec. II, we discuss the method which will be used in this paper for calculating geometry effects. In Sec. Ill, 
we calculate the exact solutions for ideal gases in two-dimensional boxes, three-dimensional tubes and boxes. By 
comparing Eqs. (2) and (3) with these exact solutions, we can see that the densities of states given by Eqs. (2) 
and (3) are very good approximations for discussing ideal gases in confined space. Moreover, in this section, we also 
consider the geometry effect on the Fermi energy. In Sec. IV, we calculate the exact solutions for ideal gases in annular 
containers. The result shows no boundary effects in such systems even the scales of the systems are very small. In Sec. 
V, we calculate the exact solutions for photon gases in confined space, and compare the thermodynamic quantities 
with the standard result of black-body radiation. In Sec. VI, we calculate an approximate result for phonon gases in 
confined space. Based on this result, we can analyze the boundary effect on the lattice specific heat. In Sec. VII, we 
analyze the geometry effects in realistic systems and compare the infiuence of geometry effects with the influence of 
fluctuations and interparticle interactions. The conclusions are summarized in Sec. VIII while some expressions of 
thermodynamic quantities are given in Appendix A. 

II. THE METHOD FOR CALCULATING GEOMETRY EFFECTS 

For an exact study of statistical mechanics in confined space, the geometry effect has to be reckoned in. To seek the 
exact solution of a grand potential, we need to find a method to perform the summation over all possible states because 
the energy spectrum is discrete in conflned space. In the following we will solve the sum by using the Euler-MacLaurin 
formula [5]: 

^ Fin) = / F{n)dn + -i^(0) - ^B,F\0) - ^B,F"\0) + ■■■, (4) 

71=0 "'O 

where Bi, arc Bernoulli Numbers: B2 = 1/6, B^ = —1/30, • • •. 

The Euler-MacLaurin formula converts a sum to a series. Generally speaking, the Euler-MacLaurin formula can 
be used as an approximate method to solve sums approximately. However, in some special cases, this series has only 
a finite number of non-zero terms, i.e., F(^'^+^'(0) = when v is greater than a certain value (note that B2V+1 = 
when > 0), and one can obtain the exact solution of the sum. 

Moreover, a more general expression for the Euler-MacLaurin formula is 

^ F{n) = / F{n)dn - - [F{m) - F{k)] + -f [F'{m) - F\k)] + -f [F"'{m) - F"'{k)] + ■■■, (5) 

which will be used to discuss the phonon gas in confined space. 

In statistical mechanics, one always replaces the summation over states by an integral in the thermodynamic limit 

^ 00. This treatment is equivalent to replacing the summation only by the integral on the right-hand side of 
Eq. (4) but neglecting the contribution from the rest terms. The information about the geometry of the system is 
contained in the summation over states since the structure of the spectrum of noninteracting particles is determined 
by the system geometry. The geometry information, however, will be lost after the replacement of the summation 
with the integral; in other words, the integral in Eq. (4) does not contain geometry information, that is to say, the 
rest terms in Eq. (4) correspond to geometry effects. 

It is natural to expect that if there is no boundary, there is no boundary effect at all. No boundary implies periodic 
boundary conditions; in other words, there will be no boundary effect if we apply periodic boundary conditions. For 
example, we will show that there is no boundary effect in the system of an ideal gas in an annular container. 

The prerequisite for applying the Euler-MacLaurin formula to solving the sum over states is that the energy spectrum 
of the system must be known. Whether one can obtain an exact solution or not depends on whether the series has a 
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finite number of non-zero terms, ff there are infinite non-zero F^^'^^^^ (fc), one obtains only an approximate solution. 
In the following we will calculate some exact and approximate solutions for ideal gases, photon gases and phonon 
gases in confined space. 



III. EXACT SOLUTIONS FOR IDEAL GASES IN CONFINED SPACE 



In this section, we first calculate some exact solutions for ideal gases in two- and three-dimensional confined space. 
Based on these exact solutions we then discuss the validity of the approximate method provided in Ref. [4], which 
can be used to calculate the geometry effect in the system with an irregular boundary. Moreover, as an application, 
we discuss the boundary effect on the Fermi energy. 



1. Exact solutions for ideal gases in two-dimensional cases 

In this part, by using the Euler-MacLaurin formula, we calculate the boundary effects on ideal gases in two two- 
dimensional cases exactly: a) ideal gases in two-dimensional boxes, and b) ideal gases on the lateral surface of 
cylinders. 



a) Exact solutions for ideal gases in two-dimensional boxes 



Next we discuss the exact solution for ideal gases in two-dimensional boxes by using the Euler-MacLaurin formula. 
The energy spectrum of a particle in a two-dimensional rectangular box of sides and Ly is 



E{nx,ny) 



Then the grand potential of the system can be expressed as 

InS = T £ £ In [l T ze-f^^^"^'"^^' 



ria- — 1 riu — l 



(6) 



(7) 



In this equation and following, the upper sign stands for bosons and the lower sign for fermions. By use of the 
Euler-MacLaurin formula, we can perform the summation exactly: 



In 



OO OC 



dnxduy In 



1 =F ze 



-(3E(n^,ny) 



±- 





OO 

1 



OO 

11' 



± ^ / drix In 





1 =F ze 



-/3B(nx,0) 



duy In 



1 =F ze' 



-l3E{0,ny 



T-ln(lT^). 







We have 



A2' 



InS = ^h2iz) - ^jh3/2{z) + ^hi{z) 



(8) 



(9) 



where S = L^Ly is the area and L = 2(L^ -|- Ly) the perimeter of the box, A = h/V^wrnkT is the thermal wavelength, 
and the function 



1 f x'^~ 



dx 



(10) 



equals the Bose-Einstein integral ga{z) or the Fermi-Dirac integral fa{z) in boson or fermion case, respectively. 

Note that Eq. (9) is the exact solution for the grand potential. The first term is just the result obtained in the 
thermodynamic limit S ^ oo; the boundary effect is described by the second term which is proportional to the 
perimeter L. Comparing Eq. (9) with the corresponding approximate result given in Ref. [4] 
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InS = ^h2{z) - \jhs,2{z) + ^hr{z), (11) 

we can see that the first two terms are completely equal to each other and the error introduced by the approximation 
appears only in the third term which describes the topology effect. The contribution from the topology effect is 
proportional to the Euler-Poincare characteristic number x, in this case x = 1 — J' = 1 since r, the number of holes in 
this two-dimensional box, equals zero. Prom this, we can see the close analogy between the two results expressed in 
Eqs. (9) and (11). This close correspondence, notice that Eq. (9) is an exact solution, can be regarded as an evidence 
of the validity of the method given in Ref. [4] . 



b) Exact solutions for ideal gases on the lateral surface of cylinders 



Consider an ideal gas confined on the lateral surface of a cylinder with radius R and length L. Such a case is just 
the same as the ideal gas confined in a two-dimensional box of sides = 2i:R and Ly = L, but the wavevector along 
X is determined by periodic boundary conditions: 



fcx = ria; = 0,±1,±2,- • •, 
while the wavevector along y is determined by fixed-end boundary conditions: 



riyTT 

Ky = , Hy = i, Z, • • • . 

The energy spectrum of a particle is 

E{n.,ny) = — {kl + kl). 
Then the grand potential can be calculated exactly by use of the Eulcr-MacLaurin formula: 

In S = T ^ 5] In [l T ^e-'5^("-"^)" 

fix ny 

= ^h2{z)-^jh3/2{z), 

where S = 2itRL is the area of the surface, and L = 2 • 27ri? = AnR is the total length of the sides. 

Furthermore, we can calculate the corresponding thermodynamic quantities directly, e.g., the specific heat 



cv 

Nk 



'hi{z) ^ho{z) 



1 



3 /i3/2(^) 1 h\^'^iz)hi/2iz) 



^QhY'iz) 8 



ho{z) 



where 



<T = 



'1 + 



1 L2 hj^^iz) 1 L hi/2{z) 



QAN S h,{z) 8^^h\'\z) 



(12) 



(13) 



(14) 



(15) 



(16) 



1 



{z)h_i/2{z) 



The following relation is used in the calculation: 

dz z hi{z) 



dT Tho{z) 



7- (17) 
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2. Exact solutions for ideal gases in three-dimensional cases 



We next calculate the boundary effects on ideal gases in two three-dimensional cases exactly: a) ideal gases in 
three-dimensional tubes with rectangular transverse cross-sections, and b) ideal gases in three-dimensional boxes. 

a) Ideal gases in three-dimensional tubes 

Consider an ideal gas enclosed in a three-dimensional tube with a rectangular transverse cross-section of sides 
and Ly. The z-component of the momentum pz is continuous since the length of the tube Lz is made sufficiently 
large. In this case, the energy spectrum of a particle in the tube can be written as 



E{nx,ny,pz) = 



2m 



Pz 

2m ' 



The grand potential is 



-OO „ _1 „ _1 



(18) 



(19) 



Here, wc have converted the summation over pz into an integral since Pz is continuous. Using the Euler-MacLaurin 
formula, we can perform the summations over Ux and Uy exactly: 



InS = -^/i5/2(^) - ^^h'iiz) + \^hs/2{z). 



4 A2 



4 A 



(20) 



Here the volume V — L^LyLz and the area of the lateral surface of the tube S — LzL where L = 2{Lx + Ly) is the 
perimeter of the transverse cross-section. 

The first two terms of the exact solution Eq. (20) are consistent with the approximate result given in Ref. [4] : 



InS = ^/i5/2(^) - ]-^h2iz) + |%/i3/2(^), 



4A2 



6 A 



(21) 



where x is the Eulcr-Poincarc characteristic number of the transverse cross-section of the tube and in this case x = 1- 
Comparing Eq. (20) with Eq. (21) we can see that, just as the two-dimensional case, the error only appears in the 
term corresponding to the topology effect. 



b) Ideal gases in three-dimensional boxes 



The energy spectrum of a particle in a three-dimensional box of sides L^, Ly and Lz is 

^2*2 ^2 



E{nx,ny,nz) 



2m [Ll^Ll^Ll ' 



"^xi '^yi fT'z — 1, 2 • • • . 



The grand potential of an ideal gas in this box is 

OO OO OO 

i-- = tE E Ei-[iT 

By use of the Euler-MacLaurin formula, we can calculate the sum exactly: 



ze 



OO OO OO 



lnS = T 



drixdnyduz In 



± - I I duyduz In 



(22) 



(23) 





OO OO 







±2 / / duxduz In 




1 =F ze' 



-l3E{nx,0,n^ 







4/ 


J duxdriy In 



1 zp 2:6"'''^^°'"'""^^ 
1 ip ^e'^^^"*'""'"^ 
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drix In 



1 ^ ze-^E{n.fi,0) 



duy In 



1 =F ze 



-/3E(0,n^,0) 



1 

^4 



duz In 



1 =F 26' 



± -ln(lT^)- 



(24) 



Then we have 



InS = ^h5/2{z) - \-^h2iz) + Y^jf^3/2{z) - ^hi{z), 



(25) 



where V = L^LyL^ is the volume, S = 2{LxLy + LyL^ + L^L^) the area of the surface and L = 4 (L^ + Ly + L^) the 
total length of the sides of the box. 

We know that a box will change to a tube if its one side, e.g., L^, is made very long. Assuming Lz S> L^, Ly, we 
then can take the approximation S ~ 2(Ly + Lx)Lz and L ~ 4^^. It is easy to see that in this case Eq. (25) is the 
same as Eq. (20) except the last term which is negligibly small. In above discussions, we only concentrate on the case 
that the temperature of the system is higher than the critical temperature. 



3. The validity of the approximate method for calculating geometry effects with irregular boundaries 

In Ref. [4], for calculating geometry effects, we first expand the grand potential InS = =F^ln(l =F ze~^'^=) as a 

series of ze~^'^' , and then sum over all terms of the expansion after calculating these terms by use of the result given 
by Kac [6] . Clearly, the validity of this treatment depends on the value of the expansion parameter ze~^'^' . Strictly 
speaking, however, when the fugacity z > 1, such as in Fermi-Dirac statistics, the grand potential InS can not be 
expanded in such a way. In the preceding paper [4] we perform the summation of the expansion of the grand potential 
without proving the convergence of the series, but, instead, we analyze the validity of the result by comparing it with 
the grand potential of an ideal quantum gas in free space and show that the grand potential in free space, whose 
validity is fully accepted, is just the zeroth-order approximation of our result. With the help of the exact solutions 
obtained above, we can make further improvement on the analysis of the validity of this treatment. 

Comparing the results obtained by the approximate method developed in Ref. [4] with the exact solutions calculated 
in this section, we can see that this approximate method is valid for calculating the geometry effects in two- and three- 
dimensional ideal gas systems. So far we have applied two methods to performing the summation in the grand potential 
for ideal gases in confined space: One is an approximate method [4] based on a mathematical work given by Kac [6] . 
Another is an exact one based on the Euler-MacLaurin formula. The approximate method applies equally well to the 
exact method in calculating the boundary effect, and the error introduced by the approximation appears only in the 
topology terms and is indeed small. 

4. The boundary effect on the Fermi energy 

Using the result obtained above, we can consider the boundary effect on the Fermi energy. Prom Eq. (25), we can 
obtain 



N = g 



(26) 



where 3 is a weight factor that arises from the internal structure of the particles (the number of internal degrees of 
freedom). 

The Fermi energy /xq is the energy of the topmost filled level in the ground state of the N electron system, so 



N 



1 (f?\"' 3/2 1,1 , 1 1/2 1 



^e^^fej ^° -4^4;^l2;;^J '^°+16^;^l2;;^) 



(27) 



The Fermi energy can be exactly calculated from this equation; however, for clarity we only calculate the next-to- 
leading-order correction: 
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= iil"^"' + (29) 

Note that, for an electron gas the spin weight is g = 2. Here /ig'^^'^ is just the Fermi energy in free space, and ^j^^°'"^'^°-^y 
reflects the influence of the boundary. It is easy to see that the existence of a boundary enhances the Fermi energy, 
and the Fermi energy increases with the decreasing size of the system. This result is not difficult to understand: In a 
finite size system, ground state energy and energy level spacing grow larger as the size of the system decreases. 

IV. EXACT SOLUTIONS FOR IDEAL GASES IN ANNULAR CONTAINERS 

The exact solutions for ideal gases in annular containers can also be calculated with the help of the Euler-MacLaurin 
formula. In this section, we consider an ideal quantum gas confined in an one-dimensional ring of radius R. The 
energy of the particle in such an annular container is determined by periodic boundary conditions: 

r^^ n = 0,±l,±2,---. (30) 



2mii2 



The grand potential of the ideal gas is 



lnS = =F J2 In (1 T ^e-'^^") 

n— — oc 

= ^ In (1 T ^) T 2 In [l T ^e-'^"''''/(2-«')] . (31) 

n=l 

The summation over the discrete parameter n can be converted to an integral exactly, leading to 

InS = t2 / rfnln 1 T ^e-''"''^'/(^™^') 
Jo L -I 



= jh/,{z), (32) 



where L = 2ttR is the perimeter of the container. There is no boundary effect because, for a ring, there is no boundary. 
This is an interest result. Generally speaking, for a finite perimeter L, the energy levels e„ are always discrete, so the 

•oo 

summation in Eq. (31) has to be performed approximately by replacing it by an integral: J2 ~^ /-oo '^^^^ 

n= — oo 

replacement is based on the assumption that the momentum is continuous, which is valid only when L ^ oo. In 
many cases, such a treatment is an approximation that neglects the influence of the boundary; however, in an annular 
container this treatment gives an exact solution. 

This result implies that, in a ring, e.g., an ideal electron gas in a conductor ring, even the perimeter is very small, 
then the energy level spacings are very large, we can replace the summation by an integral safely. The thermodynamic 
behavior of an ideal gas in a ring is the same as that in infinite free space. Generally, it is easy to prove that there is 
no boundary effect when one applies periodic boundary conditions, e.g., an ideal gas confined on a torus. 



V. EXACT SOLUTIONS FOR PHOTON GASES IN THREE-DIMENSIONAL CAVITIES 

In this section we calculate the exact solution for black-body radiation in confined space. From the electromagnetic 
field theory, it is found that the allowed values of k (the momentum p = ?ik ) in a rectangular cavity of sides L^, Ly 
and Lz are 

7r 7r TT 

kx — — Tlx^ ky = — TJy , kz = — flz, Tlx, f^yi f^z — 0) 1) 2 • • • , ('^'^) 
Lix J-iy Liz 

and the energy of a photon is 
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nt nr. 



E{nx,ny,nz) = cp = ttHcJj^ + ^ + J^^ 



(34) 



where c is the velocity of ligiit. Note that, {rix, Uy, riz) 7^ (1, 0, 0), (0, 1, 0), (0, 0, 1), or (0, 0, 0), because if any two of 
the integers nx,ny, are zero, then all of the components of the electromagnetic field will be zero. Since the number 
of photons is indefinite, then the chemical potential of the system is identically zero, we have 



InS = 2 



00 00 CX) 



na;=l ny=l nz=l 



(35) 



where 



■ln{l - exp[-/3£'(nx,nj/,n2)]} 



(36) 



and the factor 2 comes from the two possible polarizations. The summation in Eq. (35) can be performed by using 
the Euler-MacLaurin formula: 



lnH = 2 



^ "^^'^ ^(4) - -AW,m + 4t7^C(2) + Ail, 1, 0) + Ail, 0, 1) + A{0, 1, 1) 



A3r(5/2) 



4 A2 r(2) 



16 A r(3/2) 



(37) 



where, as defined in Sec. Ill, S is the area of the surface and L the total length of the sides of the cavity, C('^) is the 
Riemann zeta function, A = 2y/7Thc/{kT) is the mean wavelength of the photon. 
For the total energy in the cavity, from Eq. (37) we obtain 



(38) 



where 



Binx,ny,nz) = 



exp[(3E{nx,ny,nz)] - l" 



(39) 



It is easy to see that the first term of Eq. (38) corresponds to the Stefan law of black-body radiation: U ocT^ [7]. It 
follows that the specific heat is 



Cv = k<V 



15/i^c3 



ikTf - S 



3C(3) 



ikTf 



24fic 



ikT) + 2 [£)(!, 1, 0) + D{1, 0, 1) + D(0, 1, 1) 



(40) 



where 



Dinx,ny,n^:) 



[l3Einx,ny, n^)]^ exp[/3^(na,, n^, n^)] 
{exp[^i;(nx, riy, n^)] - l}^ 



(41) 



The contribution from the first term, which is proportional to the volume of the system, is just the result obtained 
under the hypothesis that ^ 00 (the thermodynamic limit). The last three terms, which correspond to the 
contributions from certain states, are so small that can be neglected. 

Furthermore, from the result obtained above we can learn that the density of states of a photon gas in confined 
space can be approximately written as 



p{w)(L. 



-K^C^ 4 TTC^ 



(42) 



where the second term describes the boundary effect. In other words, one can replace the summation over states 
by the integral J p{u})(L). 
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VI. PHONON GASES IN CONFINED SPACE 



In this section, we discuss the boundary effect on phonon gases. SoUds are crystal lattices of atoms, each atom 
behaving as an coupled harmonic oscillator in the harmonic approximation. Phonons correspond to the normal 
modes of the system and behave as independent quantum oscillators. Assume that the system is a rectangular solid of 
dimensions L^-, Ly, L^. The wavevector of the phonon k is restricted by fixed-end boundary conditions (not periodic 
boundary conditions because the system we considered is finite) to the values 



kx — 

ky = 

kz — 



27r 

Lx 

Lz Lz 



NxTT 



Ly 

NzTT 



(43) 



where Ni, i = x,y,z, is the number of atoms along i-axis, and the total number of atoms N = N^NyNz- Strictly 
speaking, there arc only A^j — 1, not iVi, allowed values of fc^, because the solution for Niw/Li permits no motion of 
any atom [8]; however, such an influence is so small that can be neglected safely. 
The grand potential is 



InS = -13 Uo + E 2^^0 " ^ ~ ^'^"^^'^ 



N 



N 



^ln(l - 6-'^'''=''=') - 2^1n(l - e-^"''"') 



(44) 



where 



LJi = ci \k\ , uJt = ct \k\ . 



(45) 



There are three polarizations for each value of k: two of these are transverse and one longitudinal, ct is the velocity 
of transverse waves and q the velocity of longitudinal waves. The summations in In S can be converted to a set of 
integrals by using the Euler-MacLaurin formula Eq. (5). Neglecting higher orders, we have 



InS = -■ 



+ 



kT 
2 16 



-N 



a 



{Oi + 20t) 



9 6*; 
ST 
1 



^ + 31n(l-e-*)-D3(| 
ln(l-e-*) 



(|)..,„(i-.-*)-z,.(|) 



(46) 



where 



k ' ka 



irhct 
ka 



(47) 



and a is the lattice constant. Like the Debye frequency in the Debye approximation we introduce the cutoff frequency 
LOp, which will return to the Debye frequency so that 6p will return to the Debye temperature 6d when one neglects 



the contribution from the boundary effect. As an approximation, we can take Up : 

Daix) by 



* £ 



ujd- Here we define the function 



(48) 



the case oiv = 'i gives the Debye function. 

From Eq. (46), we can calculate the internal energy 



U 



16 



1 Gi Ji 
- — + 2— 

3 It t 



(49) 



and the specific heat 
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In the high temperature Umit {T ^ 9p,6i,9t or a; <C 1), the internal energy and the specific heat are 



U : 



mkT 



[Op 
\0o 



20 V T 



16 



1 ef + 2el 



72 T2 



(51) 



and 



Cv = 3Nk / 



l-4( ^ 



20 V T 



_ 37r 
16 



1 - 



1 ef + 26>f2 

72 T2 



(52) 



where 



= 18.^^ n 



V \k 



cf c? 



(53) 



and ^0 equals the Debye temperature. 

Correspondingly, in the low temperature region {T <^ 6p,6i,9t ov x ^ 1) the internal energy and the specific heat 
become 



U = <po + Nk 



9 fe. 



rp4 



^(^.+2^,) + JC(3) 



yd 



(54) 



and 



Cv = Arfc^^--fc_C(3) 



(55) 



Each expression of the thermodynamic quantities has been divided into two parts: the first part, though it is also 
influenced by the boundary, is almost the result given by Debye; the second part, which is proportional to the area of 
the surface, describes the boundary effect. 

Moreover, similar to the case of ideal gases in annular containers, we can easily prove that there is no contribution 
coming from the term which is proportional to the area of the surface of the system if we apply periodic boundary 
conditions. 



VII. GEOMETRY EFFECTS IN REALISTIC SYSTEMS 



In this section, we shall briefly discuss the geometry effects in realistic systems. Such effects are usually nc;gligible 

in the thermodynamic limit {V > oo). The subject which will be of most interest to us will be that of boundary 

effects on the thermodynamic property of the system in which the thermal wavelength of the particle is comparable 
to the system size, e.g., mesoscopic systems or optical microcavities. 

1 ) Small size system,s: In statistical physics, one often uses the thermodynamic limit {V — > oo) as a convenient 
mathematical device in macroscopic limit. The system approaches the macroscopic limit once its size is much larger 
than the thermal wavelength. The size of a mesoscopic system is between microscopic and macroscopic scale. In 
mesoscopic systems, of course, many of rules in macroscopic physics, such as the thermodynamic limit (which neglects 
the boundary effect), may not hold. Many novel phenomena exist that are intrinsic to mesoscopic systems, e.g., size 
effects [9]. 

In realistic systems, the magnitude of boundary effects is determined by the ratio X/L, where A is the thermal 
wavelength and L denotes the linear size of the system. The thermal wavelength depends on two factors: the 
temperature of the system T and the mass of the particle. For electrons, from room temperature to lOK the thermal 
wavelength is in the range of about Anm to 2Anm. This means that in a nanosystem the boundary effect becomes 
important and can not be neglected. Furthermore, in ultra-low-temperature physics, obviously, the boundary effect 
will become remarkable. Up to now, the lowest temperature can be obtained in experiments is of the order of nK [10]. 
In this ultra-low temperature scale, the wavelength of a electron is of the order of 1mm, the wavelength of a hydrogen 
atom is of the order of 0.1mm; in other words, in such a case, the boundary effect can be measured in macroscale. 
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In the case of ideal gases confined in annular containers or crystal lattices with periodic boundary conditions, 
however, we have shown that every system, very large or very small, will not show any boundary effects, that is to 
say even in the small systems like torus C360 and C240 molecules, if we can still treat such systems by the statistical 
method, there are no boundary effects so long as the systems have the toroidal structure. Recently, there are many 
studies on nanorings, especially the metal nanorings [11]. It might be expected that the method and the result of ideal 
gases in annular containers can be used to treat the electron gases in metal nanorings. Moreover, metal wires having 
diameters in the range of nanometer are very important for nanoelectronics and other nanodevice applications, and, 
recently, ultra-thin metal nanowires have aroused growing interest in condensed matter physics [12]. The analysis 
presented above applies equally to such systems. 

Of special interest is the electron gas on a carbon nanotube. It is known that the electronic structure of a carbon 
nanotube can be either metallic or semiconducting, depending on its diameter and chirality [13]. For the metallic 
case we can treat the electrons on a nanotube as a free electron gas on the lateral surface of a cylinder. In such a 
system, along the perimeter of the transverse cross-section, even the size is very small so that the energy spectrum of 
the electrons is discrete, there is no boundary effect (sec Sec. IV); otherwise, along the nanotube axis the boundary 
effect correction has to be reckoned in when the length of the nanotube is very short. With the help of the method 
provided in Sec. Ill we can calculate the thermodynamic properties of such systems directly, e.g. the specific heat in 
Eq. (16). 

2) Long and narrow systems: Simple analysis reveals that the effect of boundary on the observable is proportional 
to L/^/S, where S is the area and L the perimeter of the system. The magnitude of the factor L/\/S is determined 
by the geometry of the system. If the shape of the system is long and narrow, in principle, the factor L/^fS can take 
on an arbitrarily large value and hence the boundary effect will become significant even in macroscale. 

3) Black-body radiation in small cavities: Prom visible light to infrared, the range of the wavelength is from 0.4 
to about 10^ micron. This means that if the size of a system is of the order of micron, the boimdary effect must 
be reckoned in. In Sec. V we have calculated the exact solution for black-body radiation in a three-dimensional 
rectangular cavity, from which we can estimate the influence of the boundary effect on black-body radiation. 

From Eq. (40) wc can estimate the influence of the boundary effect on the specific heat of a photon gas confined 
in a small cavity by comparing the contribution from the surface term (the second term of Eq. (40)) with the result 
which is obtained based on the thermodynamic limit (the first term of Eq. (40)). For this purpose, we introduce a 
ratio r] = Cy/Cy, where Cy is defined as the second term and Cy the first term of Eq. (40). For the linear size 
of the system L = 1mm and 0.1mm, when the temperature T = SOOif, the ratios rj are about 1% and 10%; when 
T = lOOK, the ratios r] are about 3% and 30%; moreover, when T = lOK, L = 1cm, the ratio r] ~ 3%. This result 
shows that, relative to the case of particles with non-zero rest masses (whose wavelengths are relatively short), the 
boundary effect is significant in photon systems. 

4) Systems consisting of sm,all grains: In Sec. VI we have discussed the boundary correction to the Dcbye theory 
of the specific heat of crystals. From Eq. (55) we can learn that at low temperature the boundary effect will become 
significant since the contribution from the boundary term is proportional to T^, whereas the term corresponding to 
the Dcbye's result is proportional to . Therefore, the contribution from the boundary effect will become more and 
more important with a drop in temperature. Of course, the smaller the system size, the stronger is the influence of 
the boundary. Consider a macroscopic system consisting of small grains. Such a granular material can be obtained, 
for example, by grinding a piece of crystal to a fine powder. In such a system the volume is almost the same as the 
volume before grinding, whereas the total area of surface, which equals the sum of the area of all grains in the system, 
becomes much larger. We can rewrite the expression of the specific heat in the following form for clarity: 

Cv = ^A{T) - ^B{T), (56) 

where the area of surface S ~ V"^^^, A{T) and B{T) are certain functions of the temperature. If we divide the volume 
V up into n cells of volume V/n, the total volume of the system is roughly speaking still V; however, the total area 

of surface becomes n(V/n)'^^^ = ^pnV'^^'^. This means that the boundary effect will become important in granular 
materials. It can be found that before grinding the correction to the standard result (the term which is proportional to 
V) is on the order a/ \jv , here we have not reckoned in the influence of temperature though it is of equal importance, 

whereas after grinding such a correction increases to the order <fn {^j \fV^ . 

5) Geometry effects vs. fluctuations: The geometry effect may be visible for small size systems. In such systems, 
however, the fluctuations also may not be neglected since the total number of particles N is not macroscopic. A 
natural question that may be asked is: Can one distinguish geometry effects from the fluctuation background? Based 
on the results calculated above, we can estimate in what cases the geometry effects will become apparent. Next we 
consider some special cases as examples. 
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The influence of fluctuations is on the order VN/N. In metals, take copper as an example, one cnr contains about 
10^^ free electrons. The ratio / ^l^^'^ reflects the influence of the boundary on the Fermi energy. When the 

linear size of the system is 5nm, the ratio /Xq°""'^'"^^/^q'^^'^ ~ 3.3%; however, in this case the fluctuation is only about 
0.9%; when the size of the system is 30n?Ti, the ratio /iQ°""'''"^^//iQ'^^'^ ~ 0.6%, but the fluctuation is about 0.06%. 
This means that the magnitude of the fluctuation, though also can not be neglected, is sufHciently smaller than that 
of the boundary effect so that one can distinguish the boundary effect in nanoscale metal systems. 

For photon gases with linear size L = 1mm (macroscalc) and temperature T = 3001^ (room temperature), the ratio 
T] (= Cy/Cy) is about 1%, but, in this case, the fluctuation is negligibly small; for the case of L = 1cm and T = IQK 
the ratio r? ~ 3%, but the fluctuation is only about 0.1%. That is to say, the boundary effect can be measured in 
macroscopic photon systems. 

6) Geometry effects vs. particle-particle interactions: In most realistic systems, there will be interactions between 
the particles. Next, we compare the contribution from the geometry effects with the influence of the interactions 
among the particles. 

(a) Bosc case: For dilute Bose gases, the weak interactions between particles and their low density allow for an 
accurate theoretical description of the effect of particle-particle interactions based on the binary collision method. 
When T > Tc, where Tc is the critical temperature, the equation of state can be expressed as [14] 

^ = ^<75/2(^)-^5l/2(^) + ---, (57) 

where a is the scattering length; the influence of the particle-particle interactions is described in the second term. For 
geometry effects, Eq. (25) gives the equation of state for an ideal Bose gas in confined space: 

PI IS 

^ = l-5.95/2(^)-^^52(2) + ---. (58) 

From this result we can estimate the contributions from the geometry effects and the interparticle interactions: 
The leading contribution from the boundary is —{l/4:)S/{VX'^)g2{z)] the leading contribution from the interparticle 
interaction is —{2a/X'^)gy2{z). By comparing these two contributions, we can see that there exists a regime where 
the geometry effect dominates in an imperfect gas. Introduce the ratio between these two contributions 

"6 = - — ^ (59) 

where a denotes the magnitude of the scattering length. The regime that the geometry effect becomes dominate can 
be determined by the condition r]b > 1- 

Using this result, we next compare the influence of the geometry effects with the influence of the interparticle 
interactions on Bose-Einstein condensation. We will show that the geometry effects, compared with the interparticle 
interactions, are important in some recent experiments on Bose-Einstein condensation in dilute atomic gases. 

In the recent experiments on Bose-Einstein condensation, the shape of the condensate is either cigar-shaped, with a 
diameter about 15/im and length SOOfxm ~ 5mm, or approximately round with a diameter of 10 ^ 50/im [15,16]. As 
an example, the Bose gas we considered here is near the critical point where the fugacity z ^ 1 and the temperature 
T ~ Tc, and the shape of the corresponding condensate is cigar-like (or thread-like), in other words, the shape of the 
system is long and narrow. In such a system, the total area of the surface is approximately equal to the area of the 
lateral surface, so we have 

dmax~0.12 — , (60) 

a 

where dmax is the maximum size of the transverse cross-section of the system in which the geometry effect becomes 
dominate. The parameters in some recent experiments have the following orders of magnitude [16]: For ^H, the 
scattering length a ~ 0.065nm, the thermal wavelength A ^ 0.25/im,, so we have dmax ^ 0.11mm; for Li, a ^ — 1.5nm, 
A ~ 1.2/xm, so dmax ~ 0.12mm; for ^^Na, CL ~ I.Snm, A ~ 0.26/xm, so dmax ~ 2.9 x 10 "^mm; for ^'^ Rb, a ~ ^.^nm, 
A ~ 0.45/im, so dmax ~ 4.3 x 10~^mm. Prom this result we can see that, compared with the interparticle interactions, 
the geometry effects can not be ignored in the experiments on Bose-Einstein condensation of and Li. For example, 
in the experiment on Bose-Einstein condensation of atomic hydrogen the condensate is 15/zm in diameter and 5mm 
in length [16], so the ratio r]b ~ 7.5, i.e., the geometry effects are more important than the interparticle interactions 
in the experiment on Bose-Einstein condensation of atomic hydrogen. 

Strictly speaking, in the recent experiments on Bose-Einstein condensation the atomic gases are trapped in magnetic 
traps rather than boxes with hard walls. Usually, such traps are approximately described as effective three-dimensional 
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harmonic wells cylindrically symmetric about the z-axis. In above analysis we assume that the Bosc gases are trapped 
in hard- wall tubes, i.e., for simplifying the treatment, we replace the harmonic wells by tubes, but the result, to orders 
of magnitude, is reliable. 

(b) Fermi case: For an imperfect Fermi gas the equation of state is [14] 

^ = {2j + 1) ^h/2{z) - 2j{2j + l)-^fi/,iz) + ■■■, (61) 

where j is the spin of the particle. Like the treatment in the Bose case, we can also compare the geometry effects with 
the interparticle interactions by comparing the second terms in Eqs. (25) and (61) directly. Next, as an example, we 

compare these two influences on the ground state of a Fermi system. 

The ground-state energy for an imperfect Fermi gas is given by {j = 1/2) [14,17] 

the contribution from the interparticle interactions described in the second term. The corresponding chemical potential 
is 

I n \ , I A 



where iiq is equal to the Fermi energy of the ideal Fermi gas plus the second term which describes the contribution 
from the interparticle interactions. Similarly, as expressed in Eq. (28), the Fermi energy containing the contribution 

r 1 1 1 J -1 . free , boundary t _i_ ji • >• boundary / interpartiele 

from tlic boundary can be written as fiQ ~ fi^ + Mo • Introducing a ratio rjf = /Mo ' "^^^ 

compare the corrections due to the geometry effects with the corrections that emerge from the interactions among 
the particles. Direct calculation gives 

Taking a box of side-length d as an example, we have ry/ = (3/8)(37r^)^/^Z^/ (ad), where I = {V/N)^^^ is the mean 
interparticle distance. Obviously, if ijf > 1, compared with the interparticle interactions, the geometry effects will 
dominate in the system. ??/ > 1 gives the condition 

J < 1-160^. (65) 

This condition is equivalent to the requirement of low densities and small system sizes. That is to say, in a small 
dilute Fermi gas, the influence of geometry effects will exceed the influence of interparticle interactions. 

Furthermore, in this paper we also consider the geometry effects on photon gases and on phonon gases. The photon 
gas can be regarded as a genuine ideal gas in statistical mechanics since the cross section of photon-photon scattering, 
which is only a loop correction, is exceedingly small. Therefore, there is no contribution from the interparticle 
interaction. In the harmonic approximation, the phonon gas can also be regarded as an ideal gas. That is to say, 
in photon gases and in phonon gases the contributions from the geometry effects are the main corrections to the 
standard results. 



VIII. DISCUSSIONS AND CONCLUSIONS 



The key problem in statistical mechanics is to calculate the partition function, i.e., to solve the sum over all 
possible states. Often, this sum is difficult to calculate. To calculate the partition function approximately, in the 
thermodynamic limit V ^ oo, one can convert the sum to an integral by introducing a state density. In so doing, 
however, the information about the system geometry has been lost because the approximation that the spectrum 
is continuous has been taken during this process. The information about the system geometry is embodied in the 
structure of the spectrum, so if we want to involve the geometry information in the thermodynamic quantities, we have 
to solve the sum over the states directly. We have developed two methods to solve the sum: One is an approximate 
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method which can be used to deal with the system with an irregular boundary [4] . Another method considered in the 
present paper is based on the Euler-MacLaurin formula; sometimes this method leads to exact solutions. However, 
the prerequisite for applying the Euler-MacLaurin formula is that the spectrum is already known. 

The method for performing the summation by the Euler-MacLaurin formula is easy to apply to the cases of ideal 
gases in external potentials and nonideal gases if the spectrum of the system is already known. Fortunately, the 
spectrimi can be obtained by using many systematic methods exactly or approximately. Such cases we will discuss in 
detail elsewhere [18]. 

Intrinsically the boundary effect arises from the interaction between the particles and the boundary. The influence 
of this interaction can be embodied in boundary conditions. The fix-end boimdary condition, like the case of a particle 
in an infinite depth potential, reflects the interaction between the particles and the wall of the box. However, if there 
is no boundary, of course, there is no interaction between the particles and the boundary; in other words, there is no 
boundary effect. In such a case we may apply periodic boundary conditions, and then the boundary terms will vanish 
in the exact solution. 

Prom the results given above, we can see that the grand potential of a system in confined space is less than that in 
free space since the sign of the leading geometry term, e.g., the second term in Eq. (20), is negative. It means that 
the existence of a boundary tends to reduce the number of states of the system, in other words, to reduce the volume 
of the phase space — the set of all possible states. This is just because, for ideal gases, in free space the spectrum is 
continuous while in confined space the spectrum gets discrete. Thus the number of modes in confined space is less 
than that in free space. 

The boundary term which is proportional to the area of the surface is always negative (the other terms corresponding 
to geometry effects are usually negligible). A natural question that may be asked is: What will happen when the 
contribution from the boundary effect become large enough, e.g., in the case of the size of the system is very small 
or the temperature is very low, so that the grand potential is less than zero? The answer is that in such a case the 
statistical method is no longer valid because the grand potential is the logarithm of the grand partition function and 
a negative grand potential means that the number of the states is less than one. 

In conclusion, by using the Euler-MacLaurin formula, the geometry effects on the statistical mechanics of various 
systems, such as ideal Bose and Fermi gases, photon gases, phonon gases, are discussed. From the exact solutions 
we can see that the grand potential can be expressed as a sum of two parts: one is just the result obtained in the 
thermodynamic limit y — > oo; another is the contribution from the geometry effect. The latter is consisted of a 
series of terms: the first term, by comparing with the result given in Ref. [4], reflects the influence of the shape 
of the boundary; the second term, which is proportional to the Euler-Poincarc characteristic number, reflects the 
influence of the topological property of the system; the rest terms are always negligible and the number of such terms 
is finite when the problem is exactly solvable, or else it may be infinite. Moreover, we discuss the validity of the 
approximate method for calculating the effects of boundary and connectivity, presented in Ref. [4] , which is based on 
the mathematical work given by Kac [6] and can be used to deal with the ideal gas system with an irregular boundary 
in two dimensions. We also analyze the possible geometry effects in realistic systems. It is hoped that the geometry 
effects can be encountered in the regime of mesoscopic scale. For photon gases, the geometry effects may be observed 
even in macroscale systems. In imperfect gases, we compare the geometry effects with the interparticle interactions 
and point out that there exists a regime where the geometry effects dominate in interacting systems. 

We are very indebted to Dr. G. Zeitrauman for his encouragement. This work is supported in part by LiuHui fund. 



APPENDIX A: THERMODYNAMIC QUANTITIES IN THREE-DIMENSIONAL BOXES 



Following general procedures, we can obtain the thermodynamic quantities for an ideal gas in a three-dimensional 
box although the geometry terms in Eq. (25) make the derivation get tedious. 
The equation of state: 

fcT " ■y^^5/2{z) - 4^'l2(^) + Y^y/i3/2(^) - g/ll(2), 



N 



'dz 



V 



1 s 

InS = ^/i3/2(^) - 4l^'^i(^) + i^^h^/2{z) - ihoiz). 



16 x' 



1 



Internal energy: 



NkT 2/i3/2(^) ^ ' N 
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where 



A = 



N 



1_S_ 
4A2 



Specific heat: 



Cy 

Nk 



4 /7,3/2(z) 4/li/2(z) -S2 

1 1 s 



(1 + A) 



2 7VA2 



h2{z) - -hi{z) , , . j- 

4 ' hi/2{z)-S2 



3 1 L 
64iV A 



h3/2{z) - hi/2{z)- 1— — - 

hi/2{z)-d2 



where 



ISX^ . . 1 iA\ , , 1X\ , , 
S2 = -—ho[z} - ——h_i/2{z) + - — h_i{z). 



4 V 



We also have 



dz 



16 V 



3 z h3/2iz) - 5i 



dT 2T hx/2{z) - 62' 
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